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By making use of the Green function concept of quanti- 
zation of the electromagnetic field in Kramers-Kronig consis- 
tent media, a rigorous quantum mechanical derivation of the 
rate of intermolecular energy transfer in the presence of ar- 
bitrarily shaped, dispersing, and absorbing material bodies is 
given. Applications to bulk material, multi-slab planar struc- 
tures, and microspheres are studied. It is shown that when the 
two molecules are near a planar interface, then surface-guided 
waves can strongly affect the energy transfer and essentially 
modify both the (Forster) short-range R~^ dependence of the 
transfer rate and the long-range R~^ dependence, which are 
typically observed in free space. In particular, enhancement 
(inhibition) of energy transfer can be accompanied by inhi- 
bition (enhancement) of donor decay. Results for four- and 
five-layered planar structures are given and compared with ex- 
perimental results. Finally, the energy transfer between two 
molecules located at diametrically opposite positions outside 
a microsphere is briefly discussed. 

FAGS numbers: 42.50. Gt, 12.20.-m, 42.60.Da, 80.20.Rp 



I. INTRODUCTION 

Intermolecular energy transfer as a fundamental pro- 
cess in many biochemical and solid-state systems has 
been of increasing interest y. It is often distinguished 
between two cases, namely (radiationless) short-range 
transfer (also called Forster transfer ||]) and (radiative) 
long-range transfer. In the former the distance R between 
donor and acceptor is small compared with the elec- 
tronic energy transfer wavelength Xa, R/Xa^^- The 
free-space transfer rate behaves as R~^, which can be 
explained by the instantaneous (longitudinal) Coulomb 
interaction between the two molecules. In the latter the 
intermolecular distance substantially exceeds the transi- 
tion wavelength, R/Xa ^ 1- The observed R~^ depen- 
dence of the transfer rate can be regarded as being the 
result of emission and reabsorption of real (transverse) 
photons. It is worth noting that in a rigorous approach 
to the problem (e.g., within the framework of the multi- 
polar formalism of QED [|[0]) the R~^ and R^'^ distance 
dependences are limiting cases of a unified theory |5| . 

When the two molecules are near material bodies, then 
the electromagnetic field felt by them can be quite differ- 
ent from that in free space and the intermolecular energy 
transfer can change accordingly. The effect has attracted 
attention, because it offers the possibility of controlling 
the energy transfer, with regard to potential applications, 
e.g., in high-efficiency light-harvesting systems, optical 



networks, and quantum computing. Enhanced energy 
transfer between molecules randomly distributed within 
a single glycerol droplet (of about 10/xm diameter) B and 
within a polymer Fabry-Perot microcavity [[7| has been 
observed. Using monomolecular layers of donor and ac- 
ceptor molecules (separated by distances of 10 . . . 20 nm) 
in planar microstructures, the dependence of short-range 
energy transfer on the local photon mode density has 
been demonstrated g. 

Calculations of the energy transfer rate have been per- 
formed in order to include the effect of bulk material 
microspheres (nO^-H, and planar microcavities |l4,f5 
The quantum theory given in Ref. [^ is based on a mi 
croscopic model that allows for both dispersing and ab- 
sorbing bulk material. In Rcfs. ||lO| , [ll|| the classical field 
generated by a donor dipole and felt by an acceptor 
dipole in the presence of a microsphere is substituted 
into the free-space Fermi's golden rule expression. A 
strictly quantum mechanical treatment that starts from 
a mode decomposition of the electromagnetic field ac- 
cording to the Helmholtz equation of the macroscopic 
Maxwell equations is given in Refs. ||l^,|lj,|l5[ . Unfortu- 
nately, the microscopic theory developed for bulk mate- 
rial 1^ becomes quite cumbersome when boundaries are 
present, and studies based on the standard mode expan- 
sion ||l^,|lj,|l^ cannot incorporate material absorption. 

In the present paper we give a rigorous derivation of 
the rate of intermolecular energy transfer in the presence 
of arbitrarily shaped, dispersing, and absorbing material 
bodies, starting from the quantized version of the macro- 
scopic electromagnetic field. The quantization is based 
on the introduction of Langevin noise current and charge 
densities into the classical Maxwell equations, which 
can then be transferred to quantum theory, with the 
electromagnetic-field operators being expressed in terms 
of a continuous set of fundamental bosonic fields via the 
classical Green tensor (sec [|l6|,0 and references therein) . 
In particular, we show that the minimal-coupling scheme 
and the multipolar-coupling scheme yield exactly the 
same form of the rate formula. It is worth noting that the 
formalism includes material absorption and dispersion in 
a consistent way, without restriction to a particular fre- 
quency domain, and applies to an arbitrary (inhomoge- 
neous) medium configuration. 

Here, we apply the theory to bulk material, multi-slab 
planar structures, and microspheres, with special empha- 
sis on media of Drude-Lorentz type. In particular, we 
show that the energy transfer can be strongly modified, 
if the two molecules are sufficiently near an interface and 
surface-guided waves at the energy transfer wavelength 



exist. Four- and five-layered planar structures are stud- 
ied, and the results are compared with recent measure- 
ments [|| . Finally, the effect of surface-guided waves and 
whispering-gallery waves in the case of the molecules be- 
ing near a microsphere is briefly discussed. 

The paper is organized as follows. In Section |l| the 
basic-theoretical concept of electromagnetic-field quanti- 
zation is outlined and the energy transfer rate is derived. 
Section III is devoted to applications, with special em- 



phasis on multi-slab planar structures, and concluding 
remarks are made in Section IV. Some deepening calcu- 



lations are given in the Appendix. 

II. BASIC EQUATIONS 

A. The Hamiltonian 

Let us consider an ensemble of point charges, interact- 
ing with the quantized electromagnetic field in the pres- 
ence of absorbing media. The minimal-coupling Hamil- 
tonian in Coulomb gauge reads |ll7| , [l8| 

H^fd^rf dujnuji^r,u;){{r,uj) 



E 



2mr 



Pa - gaA(fQ) 



+ i d^rp{r)${r)+ d'^r p{r)ip{r), 



(1) 



where f „ is the position operator and Pa is the canonical 
momentum operator of the ath (nonrelativistic) parti- 
cle of charge Qa and mass ttIq. The first term of the 
Hamiltonian is the energy of the medium-assisted elec- 
tromagnetic field, expressed in terms of bosonic vector 
fields f (r, us) with commutation relations 

[hir,u;)jlir',u;')] ^ 5kk'5{v^v')5{oj-uj'), (2) 



[h{v,u)Jk'{v\J)\ =0. 



(3) 



The second term is the kinetic energy of the charged par- 
ticles, and the third term is their Coulomb energy, where 
the corresponding scalar potential 4>{v) is given by 



(r) = / d^r 



,3 , /5(r') 



47r£o|r-r'| 



with 



/5(r) =^qJ{v~Va) 



(4) 



(5) 



being the charge density of the particles, and Eq the 
vacuum dielectric permittivity. The last term is the 
Coulomb energy of interaction of the particles with the 
medium. 



The scalar potential (^(r) and the vector potential A(r) 
of the medium-assisted electromagnetic field are given by 



V(p{v)= dwEll(r,cj)-}-H.c., (6) 
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where 



E-L(")(r,cj) = f dh'S^<-^^\r-r')E{r',Lu), 



(7) 



(8) 



with S (r) and <5"(r) being the transverse and longitu- 
dinal dyadic (5- functions, respectively, and 

E(r,c.) = Z./A ^ /'dVV^^(;V^G(r,r',c.)f(r',c.). 

V TTEo C^ J 

(9) 

Here, G{r,r',uj) is the classical Green tensor, which 
obeys the inhomogeneous, partial differential equation 



—^e(r,Lu) - V X Vx 



G(r,r',c.) = -<5(r-r') (10) 



together with the boundary condition at infinity [S{r) is 
the dyadic 5- function], with e(r, uj) = eR(r, w) -I- isi{r, oj) 
being the complex, space- and frequency-dependent per- 
mittivity. 

Let us consider the case where the particles are con- 
stituents of neutral molecules (at positions rj\/) that are 
well separated from each other. The Hamiltonian (nl) can 
then be decomposed into an unperturbed part Hq and an 
interaction part _ffi,it as follows 



H = Hn+Hi 
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(11) 



Ho^fd^rf dujhu;i^ir,uj){{r,uj) + y2HM, (12) 
J Jo ^j 



Hint — 2 2-^ VmAI' + 2_^ H 



M int • 



M 



Here, 

is the Hamiltonian of the Mth molecule. 
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(13) 



(14) 



(15) 



is the Coulomb interaction energy between the A/th and 
the Af 'th molecule, and 



HMint = y ^ 



<la^ 



VaM-^V^otM) 



+ E f^) ^'(f^M) + fd'rpMivMr) (16) 

is the interaction energy between the A/th molecule 
[charge density /5jv/(r)] and the medium-assisted electro- 
magnetic field. 

In what follows we shall restrict our attention to the 
(electric-)dipole approximation, so that Eq. dla) simpli- 
fies to 



Vmm' = £o ^dM'S^VM' - rM)dM, 



where 



IM 
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(17) 



(18) 



is the dipole operator of the Mth molecule. Disregarding 
the A^ term in Eq. dl9), which does not give rise to off- 
diagonal molecular matrix elements, making use of Eqs. 
(B)"®' ^^^ applying the dipole approximation, HMint 
takes the form of 

HMint = - duj d^rp,Mir,uj)'E{r,u;)+R.c., (19) 



where 



AM(r>^) = - 



huj 



dM.H 
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S^ir-TM) 



-dM^"(r-rAf). 



B. The transfer rate 



Let us consider the resonant energy transfer between 
two molecules A and B at positions ta and r^. The 
initial (final) state \i) (|/)) describes the excited molecule 
A (B), the molecule B (A) being in the ground state, and 
the medium-assisted field in vacuum, 

|i) = |a',6)®|{0}), E, = Ea'+Eb, (21) 

|/) = |a,6')®|{0}), Ef = Ea + Ek> (22) 

(cf . ) . Note that imposing this initial condition requires 
that the time of state preparation is sufficiently short 
compared with the time of energy transfer. Using the 
Born expansion p9| ] up to the second order perturbation 
theory, the (totaljrate of energy transfer can be given by 






(23) 



where pi is the occupation probability of the state |i), 
and 
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(24) 



with 

T = Hint + Hint - 
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-Hint, s -^ +0. (25) 



Ei — Hq + is 
Applying the decomposition (|l3|), we may write 

{f\f\i)^{a,b'\f\a',b) 

^{a,b'\VABW,b) + {a,b'\f\a',b), 

where 



(26) 



{a,b'\f\a',b) = {a,b'\ 
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\a',b). (27) 



Let us first consider the Coulomb term (a, b'\VAB W , b). 
From Eq. (p7| ) it is not difficult to see that 

(a,6'|VAB|a',6) =£(7^ dfc/fc (5ll(rB - r^) daa- , (28) 
where 

Aaa'(W) = {a{b)\dA(B)W[b')). (29) 

In order to calculate (a, 5'|T|a', 5), we make use 
of Eqs. (19) and (|20|), perform the summation 



and integrations over the possible intermediate states 
(20) |a',6'}/j(s,w)|{0}) and |a,6)/j(s,w)|{0}). After some 
calculation we derive, on applying Eq. (|9|) and the re- 
lationship p|,p7l. 



ImGfc;(r,r',w) 
= d^s — ei{s,u;)Gkm{r,s,uj)Gi^{r',s,uj), (30) 
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, (31) 



^(jJa'a — huJ + is 

where 

OJa'a = {Ea' - Ea)/h = {Ef - Eb)/h = LJb'b (32) 

and 



Aa(b) (r, ^) = ^-^ (r - rA{B) ) 



^a'a{b'b) 



d%r~rA(B)) 

(33) 



[note that A^(B)(r, Wa'Q(b'6))=5(r - r^(B))]. Recalling 
that ImG(r',r, — (jj) — — ImG(r', r, w), we may rewrite 
Eq. (|l|) as 



(a, b'\T\a' , b) 






d^r' /d^r 



dw 



huja'a — Ji^ + is sign (tj) 



or, equivalently, 

(a,5'|t|a',6) = l^ 



d^r' /d^r 



dw 



2i 



[db/6AB(r', uj)G{r', r, cj) A^lr, w)daa'] 



?iWa/a — huj + is sign (w) 
[dfcb. AB(r', u;)G(r', r, u;) A^(r, cj)da, j 
f^a'a ^ ^^ ~ *s sign (w) 



(34) 



(35) 
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X |df,/bG(rB,rA,Wa'a)daa'| (5(Wa'Q - ^b'b) • (40) 

It can be proved (Appendix ^ that the use of the mul- 
tipolar Hamiltonian p7| instead of the minimal-coupling 
Hamiltonian (|^) exactly leads to the same expression of 
the energy transfer rate. 

Let us now consider the total energy transfer rate 
according to Eq. (pq), by taking into account the vi- 
bronic structure of the molecular energy levels. Restrict- 
ing our attention to the Born-Oppenheimer approxima- 
tion and neglecting the weak dependence of the electronic 
transition-dipole matrix element on the nuclear coordi- 
nates (see, e.g., |Q), we may factorize the dipole transi- 
tion matrix elements according to 



daa'(6b') — ^A(B) 1^aa'{bb'), 



(41) 



The w-integral in Eq. (|35| ) may now be evaluated by 
means of contour- integral techniques, by taking into ac- 
count that the Green tensor is a holomorphic function of 
uj in the upper complex half-plane, which asymptotically 
behaves as ||lj 

lim — G(r,r',w) = -(5(r-r'). (36) 

We therefore close the path of integration by an in- 
finitely large semicircle in the upper complex half-plane, 
\uj\^oo, and subsequently subtract the semicircle inte- 
gral. It is easily seen that only the terms in Ayi(r, w) 
and As(r, uj) [Eq. (p3|)] which are proportional to uj con- 
tribute to the integral over the semicircle. 



(a,6'|f|a',6) 



semicircle 



^0-^ 



db'b^"(rB - rA)daa' 



(37) 



It is further seen that only the first term in the curly 
bracket contributes to the integral over the closed path. 
We thus arrive at 



{a,b'\f\a',b) 



db'b^"(rB -rA)da 



uj: 



Slf [db/bG(rB,rA,Wa'a)daa'] ■ 

Substitution of the expressions (|2^) and 
Eq. (pq) yields the transition amplitude 



(38) 
into 



(a,6'|f|a',fe) 






[db'bG(rB, r^, UJa'a)daa'] ■ (39) 



Note that the first term in Eq. (Bq) and the Coulomb 
term (|2q ) exactly cancel out. We eventually combine 
Eqs. (|24|) and ( p9| ) and find that the rate of energy trans- 
fer between the chosen states \a',b) and \a,b') reads as 



where dA(B) is the purely electronic transition-dipole ma- 
trix element of the transition between the lower and the 
upper electronic state of the molecule A{B), and Vaa'(bb') 
are the overlap integrals between the vibrational quan- 
tum states in the two electronic states of the respective 
molecule. Note that the vibrational overlap integrals take 
account of both displaced and distorted energy surfaces. 
Combining Eqs. (^3|) and (^) yields 



^I'a 






\Vb'bVaa' 



x\d*gG{YB,rA,UJa'a)<^A\ S{uJa' a ~ UJb' . 

which can be rewritten as 



w= I duJw{oj)aT{uj)(jf"{uj), 



(42) 



(43) 



where 



w{uj)^^(-^) \d*sG{rB,rA,u;)dA\\ (44) 



and 



and 



CT™(w) := ^Pa' \Vaa'\'^ S{uJa'a ^ Uj) 



_abs 



(c) 



y^Pfc [wfc'fcl" S{ujb'b - uj), 

b,b' 



(45) 



(46) 



respectively, are proportional to the (single-photon) 
emission spectrum of molecule A and the (single-photon) 
absorption spectrum of molecule B in free space each [EOJ . 
Thus, the rate of energy transfer is proportional to the 
overlap of the two spectra weighted by the square of the 
absolute value of the actual Green tensor. It is worth 



mentioning that Eqs. ([40|)-(46) apply to the resonant en- 
ergy transfer between two molecules in the presence of 
an arbitrary configuration of dispersing and absorbing 
macroscopic bodies. All the relevant parameters of the 
bodies are contained in the Green tensor. Note that the 
emission (absorption) spectrum observed in this case is 
not proportional to a°J^(uj)\a'^^(u!)] in general, as it can 



be seen from a comparison of Eq. (|45| ) with Eq. (B7). 

In particular when the Green tensor slowly varies with 
frequency on a scale given by the (relevant) vibrational 
frequencies of the molecules, then 'w{uj) is also a slowly 
varying function of frequency and can (approximately) 
be taken at the electronic energy transfer frequency 
LUA {~lub) and put in front of the integral in Eq. (43), 
thus 



w ~ 'w{uja) ct, 



where 



<7= / duj<7'^\cj)<7f%Cj). 



(47) 



(48) 



In this case, the influence of matter environment on the 
(total) energy transfer rate is fully contained in {^{uja)- 
Clearly, when the two molecules are near a resonator- 
like equipment, so that the molecule can "feel" sharply- 
peaked field resonances, then wluj) cannot be assumed to 
be a slo wly v arying function of frequency in general (see 
Section pl^ ). 

It may be interesting to compare the rate of energy 
transfer with the donor decay rate. Straightforward gen- 
eralization of the well-known formula for a two-level tran- 
sition yields, on applying the Born-Oppenheimer approx- 
imation, 



Ta 



where 



TAico) 






du;fA{L^)aT{^), (49) 



[d^Im G{rA, va, uj)dA] , (50) 



and a'^{uj) is given by Eq. (p5|). Whereas the decay rate 
is determined by the imaginary part of the Green tensor 
(taken at equal positions), the transfer rate is determined 
by the full Green tensor (taken at different positions). 
Thus, decay rate and transfer rate can quite differently 
respond to a change of the environment. 



III. APPLICATIONS 

A. Bulk material 

Let us first consider the case when the two molecules 
are embedded in bulk material of arbitrary complex per- 
mittivity e{uj). Using the well-known expression of the 



bulk-material Green tensor G^"''^(r, r', uj) (see, e.g., [Q), 
application of Eq. ( [44| ) yields 



w(w) 



where 
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d^G^"'^(rB,rA,c^)dA , (51) 



d^G^""^(rs,r^,c^)d 



q{uj) 



- i d*BdA - 3 
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An 
d*gR d^R 



ex.p[iq{Lu)R] 
1 



q^Lj)R^ q^Uj)R^ 
1 



q{u;)R 



with 



qiuj) = Ve(tj)- 



R = rs - r^ . 



(52) 



(53) 



The energy transfer rate is then obtained according to 
Eq. (Es). Obviously, the Green tensor of bulk material 
can be regarded as being a slowly varying function of 
frequency, so that the approximation (B^) applies. 

From Eqs. ( |5l| ) and ( p2[ ) it is seen that the en- 
ergy transfer rate includes both the small-distance case 
(Forster transfer), with the rate being proportional to 
R~^, and the large-distance (radiative) case, where the 
rate becomes proportional to R^^. Note that the expo- 
nential I exp[i(7(a;)i?]|^ = exp[— 2wni(a')i?/c], which typi- 
cally arises from material absorption, drastically dimin- 
ishes the large-distance energy transfer [y^e{uj) = n(a;) 
= n]x{uj) + ini{uj)]. In Eq. (M) local-field corrections are 
ignored. They may be taken into account by applying, 
e.g., the scheme used in Ref. [ll^ for correcting the rate 
of spontaneous decay. 

It is worth noting that the above given result, which 
is based on the quantization of the macroscopic Maxwell 
field for given complex permittivity, exactly corresponds 
to the result obtained in Ref. [pi within the framework of 
a fully microscopic approach on the basis of some model 
medium coupled to the radiation field and a heat bath. 
Already from the study of the spontaneous decay of an 
excited atom near an interface pl| it is clear that in 
the case of inhomogeneous media (of complicated atomic 
structure) a microscopic approach would be rather in- 
volved and closed solutions would hardly be found. 



B. Multi-slab planar structures 

Let us consider a planar multi-slab structure and as- 
sume that the two molecules are in the same slab. The 
relevant Green tensor (for the energy transfer between 
the two molecules relevant) of an inhomogeneous system 
of this type can always be written in the form of 

G(rs,r^,^) = G'^"'^rB,rA,c^) + G"-<=«(rB,r^,c^), (54) 
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1. Geometry of the multi-slab planar structure prob- 



where G'^""^(rs, r^,^) is the Green tensor according to 
Eq. ([5^), with e{u;) being the permittivity of the slab 
in which the two molecules are located, and the reflec- 
tion term G''®''(rB,r^,a;) insures the correct boundary 
conditions at the surfaces of discontinuity. Clearly, a de- 
composition of the type of Eq. (|4|) is also valid for other 
than planar systems, provided that the two molecules are 
located in a region of space-independent permittivity. 

To be more specific, let the z-direction be the direction 
of variation of the permittivity of the multi-slab system 
and assume that ta and tb arc in the jth slab of thickness 
dj (Fig. 1^). The reflection term in Eq. ( |5^ ) can then be 
given by 111 (see also Ref. ||) 



G'-'=f^(rB,r^,c^) 



An 
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g./3,d,grcfi(j.^^j.^^^^. 



(55) 



[fcj == y^ej{uj)uj/c; Pj = (fc^-fcp^/^]. Choosing the coor- 
dinate system such that Ry = 0, the nonvanishing com- 
ponents of G™f* read 
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-^C^[jo(fc||i?.)-( + )J2(fc||i?.)] 

hC; [Jo(fc||i?x) +(-) J2(fc||i?.)] , (56) 
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(57) 
(58) 



[J„(x) - Bessel function], where 

J.1 ^ifij(zA+ZB-dj) ^ 
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r1 p-il3j{zA+ZB~dj) 



-)2rlr'Lcos{l3jR,)e'^^'^' 
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j,1 git3j{zA+ZB-dj 



-ifjj{zA+ZB~dj) 



')2ir\r1sin{(3jR,)e"^^'^^ 



D-,\ 



D„ 



^q „<i „2i(3jd, 



(59) 

(60) 
(61) 



Here, q=p{s) means TM(TE) polarized waves, and ^'l/-, 
are the total reflection coefficients at the upper (lower) 
stack of slabs [j' < j (f>j)] of the waves in the jth slab 
(for details, see Ref. 123]). Note that when r^ and tb are 
in the top (bottom) slab, then Eqs. (55)-(p|) (formally) 



apply provided that t"^/ -, ~ and dj — are set. 

If the frequencies of the vibronic transitions that are 
involved in the energy transfer are sufficiently far from 
a medium resonance, so that material absorption (in the 
jth slab) may be disregarded, then the permittivity may 
be considered as being real and positive. In this case, it 
may be useful to decompose the integral in Eq. (Bq) into 
two parts, 

G--^" (rs , r^ , c^) = Gr"" (rs , r A , c^) + G^2<=« (rs , r^ , c^) , 

(62) 



Gi<=«(rB,r^,c^) 
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^""ai,.|«.,G»"(„,,„,,,,,,, ,63) 

'^Pj 



Gl^''{rB,rA,u;) 
= ^r »e-l/^.KG-\r.,r.,.,.„). (64) 

Obviously, G\'^^ results from waves that have a propa- 
gating component in the z-direction, whereas the waves 
that contribute to 62'^^ are purely evanescent in the z- 
direction. 



1. Interface 

Let the two molecules be embedded in a half-space 
medium (medium 1) and assume that in the relevant 
frequency interval the permittivity of the medium ei(w) 
can be regarded as being real and positive. When the 
molecules are near the interface between the two half- 
space media such that ki(zA -I- z^) ^ 1, it can be proved 
that Eqs. (|55|)-(|5|) reduce to {kiR^ < 1) 
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FIG. 2. The electronic part of the rate of energy trans- 
fer [Eq. ( |44[ )] between two molecules near a planar dielec- 
tric half-space is shown as a function of the transition fre- 
quency for z-oriented transition dipole moments and a sin- 
gle-resonance Drude-Lorentz-type dielectric [_Ra; = 0.015 At; 
za = z_g = 0.02 At; a;p=0.5a;T; 7/a;T = 10~* (solid line), 
10"^ (dashed line), and 10~^ (dotted line)]. The inset shows 
the electronic part of the corresponding donor decay rate 
[Eq. (§)]. 
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[£2('-^), complex permittivity of medium 2]. Note that for 
r^^rs, Eqs. (65)-(|67|) just give the Green tensor whose 
imaginary part determines the influence of the interface 
on the rate of spontaneous decay of a single molecule 
[£l|,|^. (For some special cases, see also Ref. Q.) Un- 
der the assumptions made, the main contribution to G'^'^^ 
comes from ^2^. Hence surface-guided waves (including 
decaying waves) play an important role and can notice- 
ably influence the resonant energy transfer. In particular 
when medium 2 is a metal or a dielectric with £2R < 0, 
(and typically £21 ^ |£2r|), then a strong effect is ob- 
served for £2r(cj) = — £i(ijj), which is nothing but the 
condition for best excitation of surface-guided waves P5| . 
In the numerical calculation of w(w) [Eq. (Q)], which 
contains the relevant information about the influence of 
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FIG. 3. The electronic part of the rate of en- 

ergy transfer [Eq. (0); units (|dAdi3|tj^/(feoc^))^/(87r)] be- 
tween two molecules near a planar dielectric half-space 
is shown as a function of the intermolecular distance 
for z-oriented transition dipole moments and a sin- 
gle-resonance Drude-Lorentz-type dielectric [a; = 1.062ljt; 
za = 2s = 0.02 At; a;p = 0.5LJT; 7/ajT = 10~^ (solid line) 
and 10~^ (dashed line)]. The dependence of w on 
the molecule-interface distance is illustrated in the inset 
[)/uT = 10"''; ZA^ZB^ 0.02 At (curve 1), 0.03 At (curve 2), 
and 0.05 At (curve 3)]. For comparison, the free-space result 
is shown (dotted lines). 

the interface on the rate of energy transfer [see Eqs. (H3)- 
([iq)], we have assumed that the two molecules are situ- 
ated in vacuum [£1(0;) = 1] above a half-space medium 
of Drude-Lorentz type and restricted our attention to a 
single-resonance medium. 



£2(0;) = £(cj) = 1 + 



lUJ^ 



(68) 



Here, wp corresponds to the coupling constant, and 
cjt and 7 are respectively the medium oscillation fre- 
quency and the linewidth. Recall that the Drude- 
Lorentz model covers both metallic (wt = 0) and dielec- 
tric (wt 7^ 0) matter and features a band gap between 
(jJt and [jJl = -^/o;^ -I- Wp . We have performed the cal- 
culations using the exact Green tensor [Eqs. (p^)-(|6l|)]. 
Comparing the results with those obtained by using the 
approximately valid Green tensor [Eq. (p4f) together with 
Eqs. (|65|)-(|67|)], we have found good agreement. 

The behavior of w{lS) is illustrated in Fig. g. It is 
seen that outside the band gap (w < lot) where £r, > the 
modification of w(yi) due to the presence of the interface 
is small even for small distances of the molecules from 
the interface. Since in this frequency domain w(ui) may 
be regarded as being slowly varying on a frequency scale 
defined by the vibrational frequencies of the molecules, 
Eq. ( |47| ) applies. Thus, the energy transfer rate is simply 
proportional to w{u:a)- 

Inside the band gap, however, the interface can sig- 
nificantly affect w[u:) if, according to Eqs. (p5h-(p7h. 
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FIG. 4. The electronic part of the rate of energy trans- 
fer [Eq. ( |44[ )] between two molecules near a planar dielec- 
tric half-space is shown as a function of the distance of the 
molecules from the surface {za = zb) for z-oriented transition 
dipole moments and a single-resonance Drude-Lorentz-type 
dielectric [ij^, = 0.85 At; cop — 0.5 lot; ■y/iOT = W~'^]- For com- 
parison, the results that are obtained by taking into account 
in Eq. (J6|) only Gf^ (dashed line) or Gf^ (dotted line) are 
shown. 

£r((jj) ~ — 1 (cj ~ 1.06 wt in Fig. g), that is to say, if the 
energy transfer transition under consideration is tuned to 
a surface-guided wave. Note that a negative real part of 
the medium permittivity can easily be reahzed by met- 
als. Careful inspection of the contributions G™'^ and 
G'^'^^ to G reveals that the enhancement of w{uj) results 
from G™^, whereas the reduction reflects some destruc- 
tive interference of G™'^ and G'^°^. Another interest- 
ing feature is that the reduction of w(a;) can go hand in 
hand with an enhancement of the corresponding quantity 
f^(w) [Eq. (|^] for the donor decay rate Ta [Eq. ^] 
(see the inset in Fig. g). 

Further, Fig. reveals that with increasing material 
absorption (i.e., with increasing value of 7) w{u!) varies 
less rapidly inside the band-gap region, and enhance- 
ment and reduction are thus less pronounced. Clearly, 
the strong influence on w{uj) of the interface which is ob- 
served for small material absorption must not necessar- 
ily lead to a correspondingly strong change of the energy 
transfer rate, because of the integration in Eq. (Esh. Nev- 
ertheless, the results show the possibility of controlling 
the resonant energy transfer by surface-guided waves. 

Figure 1^ illustrates the dependence oi w{uj) on the in- 
termolecular distance for the case when uj corresponds to 
a surface-guided wave frequency and a noticeable change 
oi'w{oj) is observed {co = 1.062 wt in the figure). It is seen 
that the R~^ dependence, which is typical of the Forster 
transfer in free space, is observed for much shorter in- 
termolecular distances. The relative minima oi w{lu) be- 
low the free-space level, which are observed for some- 
what larger intermolecular distances, again result from 
destructive interference between G^^'^ and G'^"^. Eventu- 



ally, the large-distance reduction of w{u!) below the free- 
space level results from material absorption. As already 
mentioned, the behavior of 'w{uj) in Fig. g is dominated 
by surface-guided waves that decay exponentially along 
the ±z-directions. With increasing material absorption 
the penetration depths decrease, so that on average w{uj) 
becomes closer to the free-space level. The possibility of 
controlling the ultrashort-range energy transfer by vary- 
ing the distance of the molecule from the surface is illus- 
trated in the inset. 

In Fig. the dependence of iu{uj) (again for 
a;~ 1.062 wtT on the molecule-surface distance is plot- 
ted, and the contributions to w{lu) from ordinary waves 
having a propagating component in z-direction {Gi) 
and surface-guided waves {G2'^^) are shown. It is clearly 
seen that when the two molecules are very near the sur- 
face, then energy transfer between them is mediated by 
surface-guided waves, whereas for larger distances ordi- 
nary waves play the dominant role. Note that the oscil- 
latory behavior is typical of the latter case. Clearly, for 
very large distances {za, zb 3> At) the free-space behav- 
ior is observed. 



2. Comparison with experiments 

Recently, experiments have been carried out to 
study the transfer of excitation energy between dye 
molecules confined within planar optical microcavi- 
ties [@|. In the experiments, donors (Eu'^+ com- 
plex) and acceptors (l,r-dioctadecyl-3,3,3',3'-tetra- 
methylindodicarbocyanine) embedded within a transpar- 
ent material (22-tricosenoic acid) bounded by no (weak- 
cavity structure), one (half-cavity structure), or two (full- 
cavity structure) silver mirrors are considered. To com- 
pare the experimental results with the theoretical ones, 
we have modeled the half-cavity structure by a planar 
four-layered system and the full-cavity structure by a 
five-layered system. The former consists of vacuum, di- 
electric matter (22-tricosenoic acid, e = 2.49 [Em, thick- 
ness d), metal (silver, £ = — 16.0 -I- 0.6i [g^l? thickness 
25 nm) , and vacuum, and the latter consists of vac- 
uum, metal (silver, thickness 20 nm), dielectric matter 
(the same as above, thickness d), metal (silver, thickness 
25 nm), and vacuum. In each system, the donor is situ- 
ated in the middle of the dielectric layer, while the posi- 
tion of the acceptor is shifted towards the silver mirror of 
25 nm thickness. The Green tensors of the two systems 
can be calculated according to Eqs ([5^)-(|6l|). Assigning 
to silver a Drude-Lorentz-type permittivity [p7| , it can 
be proven that in the relevant frequency interval (of over- 
lapping donor emission and acceptor absorption spectra) 
w{uj) [Eq. (g|)] and Va{uj) [Eq^^(|ol)] sufliciently slowly 
vary with w, so that [cf. Eq. ([47|)] w~ w[ijja) and, sim- 
ilarly, Ta^ Ya{m.>a)- Thus, w{ijJa) and Ta{^a) can be 
viewed as measures of the energy transfer rate and the 
donor decay rate, respectively. 
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FIG. 5. The electronic parts of the donor decay rate 
(a) and the donor-acceptor energy transfer rate (b) (aver- 
aged over the dipole orientations) of molecules in cavity-like 
systems are shown as functions of the cavity length for the 
four-layered system (dashed line) and the fi ve-layered sys- 
tem (full line) considered in Section [IIB2 (AA=614nm; 
R=~R^ = 2Anm). 



Figure ^ shows the dependence on d of fA{<^A) 
and w{uja) (averaged over the dipole orientations). 
From Fig. ||(a) it is seen that at d/Ayi~0.21 (i.e., 
d~130nm for AA = 614nni) the ratio of the donor 
decay rates for the five- and four-layered systems is 
T A{^A)\bl^ A{^A)\'^^'^■'i^ which (within the measure- 
ment accuracy) is in sufficiently good agreement with 
experimental result (see Fig. 2D in Rcf. Q). Note that 
in the vicinity of d/\A ^0.21 the ratio of the two rates 
sensitively responds to a change of d/A^. 

Comparing TAi^JA) [Fig. |(a)] with w{uoa) [Fig. |(b)], 
we see that for the four-layered system and d/\A ~ 
0.16.. .0.33 (i.e., d~ 100...200nm for AA = 614nm) both 
TA{t^A) and w{!jOa) decrease with increasing d and an 
approximately valid linear relation between the energy 
transfer rate and the donor decay rate can be established 
in agreement with experimental results in Ref . [0 . From 
the data reported in Ref. js) it could be expected that the 
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FIG. 6. The electronic part of the donor-acceptor en- 
ergy transfer rate (averaged over the dipole orientations) of 
molecul es in th e five-layered cavity-like system considered in 
Section III B 2 is shown as a function of the cavity length for 
various values of the intermolecular distance [AA=614nm; 
7?z = — 24nm; Rx=Q (solid line), 10 nm (dashed line), and 



20 nm (dotted line)]. 

linear relation between the two rates is generally valid. 
This is of course not the case. Since the energy transfer 
rate is determined by the full (two-point) Green tensor, 
whereas the donor decay rate is only determined by the 
imaginary part of the (one-point) Green tensor, the two 
rates can behave quite differently, as it is demonstrated 
in Fig. ^ In particular, the increase of the donor decay 
rate at the cavity resonances can be accompanied with a 
decrease of the energy transfer rate, because of destruc- 
tive interferences. 

In the experiments in Ref. |@], the measurements are 
performed on an ensemble of donors and acceptors whose 
distance is fixed in the z-direction but variable in the 
x-direction (Ai?2,~lnm). The question thus arises of 
whether the measured data refer to a single nearest- 
neighboring donor-acceptor pair {Rx — 0) or not. In 
Fig. H we have plotted the dependence on d of w{uja) 
(averaged over the dipole orientations) for the five- 
layered system and various values of i?^, with R^ be- 
ing fixed. We see that the rates of energy transfer be- 
tween molecules whose distances are larger than that of 
nearest-neighboring molecules can be quite comparable 
with those of the latter. Moreover there are also cases 
where the energy transfer rate increases with the donor- 
acceptor distance. The experimentally determined en- 
ergy transfer rates are thus averaged rates, which not 
necessarily show the characteristic features of single-pair 
transfer rates. Averaging in Fig. o wiujA) over all val- 
ues of Rx , the resulting curve is expected to be substan- 
tially flatter than the solid- line curve {Rx^Q), particu- 
larly when d sweeps through \a- 

An analysis of the contributions of G^^^ [Eq. (|63|)1 and 



G"-:^^ [Eq. (H)] to C"" [Eq. (|6|)] reveals that for cavity 



lengths of (i/Ayi< 0.16 (i.e., d< lOOnm for AA = 614nni) 
evanescent waves dominate the influence of the cavity 
system on both the rate of intermolecular energy trans- 
fer and the donor decay rate and lead to a strong increase 
of them. Whereas for cavities lengths of cI/Xa^O.SI 
(i.e., d^ 500 nm for A^ = 614 nm) evanescent waves only 
weakly affect the donor decay rate, they can strongly 
affect the intermolecular energy transfer up to cavity 
lengths of a few micrometers. Note that the resonance 
lengths seen in Fig. o originate from propagating waves. 



C. Microsphere 

Microspheres have been of increasing interest, because 
of the whispering-gallery (WG) and surface-guided (SG) 
waves, which may be employed, e.g., for reducing the 
thresholds of nonlinear optical processes [ p8|p9| ]. Inter- 
molecular energy transfer in the presence of microspheres 
has been considered for molecules near a small metallic 
spheroid (spheroid's linear extension -C Xa) in the non- 
retardation limit, for molecules embedded within a di- 
electric microsphere [OJl2J, and for the case where one 
molecule is inside a dielectric microsphere and the other 
outside it |13[. Here we restrict our attention to the in- 
fluence of WG and SG waves on the energy transfer be- 
tween two molecules outside a microsphere, taking fully 
into account retardation effects. 

Let £i(w) and £2(1-0) be respectively the permittiv- 
ities outside and inside the sphere. If the transition 
dipole moments are parallel to each other and tan- 
gentially oriented with respect to the sphere, the rel- 
evant (spherical-coordinate) components of G'^'^^ are 
(0A = 0B =0,0^-0) 



_ iki ^ (2M-1) J »^M,(i) 



\l{l + 1)P/(cos6'b) - cos 6(5 P/ (cos 6*5 
]^[kirAhl '{kirA}]' 



B 



kiTA 



[kirBhf\kirB)]', 



kirs 



P/ (cos 0s) 



(69) 



(for the Green tensor of a sphere, see, e.g., pOJ), and 
for radially oriented dipoles the relevant components are 

{cl)A^(f>B^O,eA^O) 



ik^ 



GZ%ArB,rA,u) = '-^Y. 



1{1 + 1){21 + 1) 



1=1 



TATE 



X Bf hY' {kirA)hY\kirB)Pl{cost 



(70) 
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FIG. 7. The electronic part of the rate of energy transfer 
[Eq. (M)] between two molecules (at diametrically opposite 
positions) near a microsphere is shown as a function of fre- 
quency for radially oriented transition dipole moments and 
a single-resonance Drude-Lorentz-type dielectric [a = 2 At; 
rA~rB= 2.02 At; i^-p — 0.5 cjt; 7/cjt ~ 10~*]. 



Bf 



(c) 



[a2Ji(«2)] ji{ai) - [aijijai)] ji{a2) 
[a2Ji{a2)]' hf' (ai) - ji{a2)[aihf'{ai)]' 



(71) 



B^ 



(^) 



£1 



{^)ji{a2)[aiji{ai)] - £2(w)i;(ai)[a2Ji(«2)] 



:i{uj)ji{a2)[aihf\ai)]' - e2(w)[a2J/(a2)]'/ii (ai) 

(72) 



where 



[ai,2 = fci,2a; a, microsphere radius; ji^z), spherical 
Bessel function; K^ (z), spherical Hankcl function; 
Pl"'{x), associated Legendre function]. 

In Fig. |7| the dependence on frequency oi w{uj) is illus- 
trated for the case where vacuum is outside the sphere 
and the two molecules are placed at diametrically op- 
posite positions (r^=:— r^), with the transition dipole 
moments being radially oriented. It is clearly seen that 
the energy transfer can greatly be facilitated at the posi- 
tions of the sphere-assisted field resonances, the enhance- 
ment of w{(jj) at the positions of SG resonances (inside 
the band gap) being larger than those at the positions of 
WG resonances (outside the band gap). Maximum values 
of w{uj) are observed where the SG resonances overlap. 
The energy transfer rate for tangentially oriented dipoles 
(not shown) is in general smaller than that for radially 
oriented dipoles. Note that when {^{lo) is sharply peaked 
at the sphere-assisted field resonances, such that it is not 
slowly varying in the frequency interval where the (free- 
space) donor emission and acceptor absorption spectra 
overlap, then it cannot be taken at the electronic en- 
ergy transfer frequency and put in front of the integral 
in Eq. (Esh. In this case, the change of the energy transfer 
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rate will be less pronounced than it might be expected 
from the frequency response of the electronic part, be- 
cause of the frequency integration. 



IV. CONCLUSIONS 
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We have given a rigorous, strictly quantum mechanical 
derivation of the rate of intermolecular energy transfer in 
the presence of dispersing and absorbing material bod- 
ies of arbitrary shapes, showing that both the minimal- 
coupling scheme and the multipolar coupling scheme lead 
to rate formulas of exactly the same form. The depen- 
dence on the material bodies of the energy transfer rate 
is fully expressed in terms of the Green tensor of the 
macroscopic Maxwell equations for the medium-assisted 
electromagnetic field. In the macroscopic approach, the 
dispersing and absorbing material bodies are described, 
from the very beginning, in terms of a spatially varying 
permittivity, which is a complex function of frequency. 
The macroscopic approach has - similar to classical op- 
tics - the benefit of being universally valid, without the 
need of involved ab initio microscopic calculations. In so 
far as such calculations for simple model systems have 
been performed, the results agree with those obtained 
from the microscopic approach. Clearly, macroscopic 
electrodynamics is valid only to some approximately fixed 
length scale which exceeds the average interatomic dis- 
tance in the material bodies. 

Whereas the donor spontaneous decay rate is deter- 
mined by the imaginary part of the Green tensor in 
the coincidence limit, the donor-acceptor energy transfer 
rate depends on the full two-point Green tensor. Hence, 
the decay rate and the energy transfer rate can be af- 
fected by the presence of material bodies quite differently. 
Our calculations for planar multilayer structures have 
shown that enhancement (inhibition) of spontaneous de- 
cay and inhibition (enhancement) of energy transfer can 
appear simultaneously. They have further shown that 
surface-guided waves can strongly affect the energy trans- 
fer, thus being very suitable for controlling it. 

In free space it is often distinguished between two 
limiting cases, namely the short-distance nonradiative 
(Forster) energy transfer and the long-distance radiative 
energy transfer. The former is characterized by the R~^ 
distance dependence of the transfer rate, and the lat- 
ter by the R~^ dependence. In particular, in the short- 
distance limit the energy transfer rate rapidly decreases 
with increasing distance between the molecules. This 
must not necessarily be the case in the presence of mate- 
rial bodies, because of the possibly drastic change of the 
dependence on the distance of the energy transfer rate. 
So, our calculations for planar multilayer structures have 
shown that the energy transfer rate can also increase with 
the distance. 



APPENDIX A: DERIVATION OF THE 

TRANSFER RATE IN THE 
MULTIPOLAR-COUPLING SCHEME 

The multipolar-coupling Hamiltonian can be obtained 
from the minimal-coupling Hamiltonian by means of the 
Power-Zienau transformation |^,y , 



n = U^'HU, 



where 



U = exp 



M •' 



(Al) 



(A2) 



with 



PM(r) = ^ fcM (r^M - I'm) 



j dX5[ 
Jo 



r-rM-A(faj^, -rjv/)] (A3) 



being the polarization associated with the Mth molecule. 
Using H from Eq. (|l|), we derive (see, for details, ||l7|| ) 

-».= fd^r f dLuhLj{''{r,Lu){{r,Lu) 
^f^2m„,a 

pi ^ 2 

+qaM / dA A {TaM -I'm) X B [rAf + A (?«„ -tm) 







M •' 



M 



-^PM(r)PM(r) 



PM(r)E(r) 



(A4) 



where B(r) = Vx A(r) [with A(r) from Eq. (^], and 



E(r) 



dcjE(r,tj) 4-H.c. 



(A5) 



and neutral molecules with non-overlapping charge distri- 
butions are again assumed. Note that in the multipolar- 
coupling scheme the operator of the electric field strength 
is defined according to 

E{v) = -1 [A(r), Til - V(^(r) - V0(r) , (A6) 

in L J 
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which imphes the following relation between E(r) and 
i{r): 



£oE(r) = £o^(r) + ^PM(r). 



(A7) 



M 



Hence, £oE(r) has the meaning of the displacement field 
with respect t o th e molecular polarization. 

From Eq. (A4) it is seen that the molecules now 
interact only via the medium-assisted electromagnetic 
field. In par ticular, in the (electric-)dipole approxima- 
tion Eq. (A4) simplifies to 



Ti. — Tio + Hi, 



(A8) 



where 



Ho^ [d^r f dujhuji''ir,uj){{r,uj) + y2'^M (A9) 
■J Jo M 



with 



^^^ = E^pL+ /d3r-^PM(r)PM(r) (AlO) 

is the unperturbed Hamiltonian of the medium- assisted 
electromagnetic field and the molecules, and 



'lint — 



M 



7iLf^=-^dME(rM) (All) 



M 



is the interaction energy between them. 

Comparing the multipolar-coupling energy given by 
Eq. (l All I ) with the minimal-coupling energy Hint given 
by Eq. (|lj) together with Eqs. Q-®, we see that the 
two energies (formally) become equal to each other, if 
we remove in the latter the Coulomb term and replace 
— [cIm, Hm]/^-^ with (Im- Having these changes in mind, 
we now fo llow step by step the derivation of Eq. ( pj| ) 
in Section [IB . Starting from the corresponding eigen- 



states of the unperturbed multipolar-coupling Hamil- 
tonian (instead of those of the unperturbed minimal- 
coupling Hamiltonian), it is not difficult to see that the 
result is again Eq. (|o|). It should be pointed out that 

the above mentioned difference between E(r) and £{r) 
[Eq. (A7)] does not affect the energy transfer rate. 



APPENDIX B: SINGLE-MOLECULE EMISSION 
SPECTRUM 

In the electric-dipole approximation and the rotating- 
wave approximation, the Hamiltonian for a single 
molecule (at position ta) that (with regard to the vi- 
bronic transitions |a') <-> \a)) resonantly interacts with 
the medium-assisted electromagnetic field reads, by ap- 
propriately specifying Eqs. (|l|)-(|20), Q 



H — Ho + Hint , 



(Bl) 



Ho= f d^r f duj huj f ^ (r, w)f (r, u) 

+ ^?iWa|a)(a| +y^huja'\a'){a'\, 



(B2) 



Hu 



a,a' 



[|a')(a|E(+)(rA)da,, +H.C.], (B3) 



where E*^"*") (r) is the positive-frequency part of E(r) de- 
fined by Eq. (A_5), and the vibronic transition-dipole ma- 



trix elements da'a of the vibronic transitions are given, 
in the Born-Oppenheimer approximation, by Eq. (|4l|). 
Let us assume that the molecule is initially (at time t — 
0) prepared in a statistical mixture of vibrational states 
in the upper electronic state and the medium-assisted 
electromagnetic field is in the vacuum state, i.e.. 



p(i = 0) = 5]p,n«'>(«'l®l{0})({0}| 



(B4) 



The time-dependent spectrum of light observed at posi- 
tion r (in free space) by means of a spectral apparatus of 
sufficiently small passband width can be given by (see, 
e.g., 11) 



5(r,ws,r)= / di2 / dii 



^-iuJs{t-2-ti) 



^0 

x(E(-)(r,t2)E(+)(r,ii))], (B5) 

where loq and T are respectively the setting frequency 
and the operating time of the spectral apparatus. In 
order to calculate the electric-field correlation function 
associated with the light emitted by the molecule during 
the spontaneous decay of the upper electronic state, we 
may restrict our attention to the perturbative expansion 
of the time evolution operator up to the first order in 

iJint 111, 



ih 



(B6) 



We make use of Eqs. (B3), ( A5) [together with Eq, 



(B4), and (B6), apply Eq. (B5), and derive after some 
calculation, on recalling the relation (30), (see also [p2|) 



lim T-^S{r,ujs,T) 

T — >oo 



2TrJ2Pa'\Va'a\'^\Fir,rA,LUa'a)fS{uJs~UJa'a), (B7) 



where 
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F{r,r A, UJa'a) 

1 f°^ up' 
= / duj^yImG{r,rA,uj)dAC{oJa'a~(^) 

TTEo Jo C 

^■G{r,rA,uja'a)dA (B8) 



- 2 



[C(x) = nSjx ) + iV/x; V, principal value]. In the deriva- 
tion of Eq. (B7), retardation has been ignored and the 
relation 

lim - [ dt2 f dti e-"^^*^-*^^ 
T^oo T Jo Jo 

,2/ 



sinVr/2) ^ 



(B9) 



has been used. 



* on leave from the Institute of Physics, National Center for 

Sciences and Technology, 1 Mac Dinh Chi Street, District 

1, Ho Chi Minh city, Vietnam. 
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